Generation of multi-photon entanglement by propagation and detection 
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We investigate the change of entanglement of photons due to propagation. We find that post- 
selected entanglement in general varies by propagation and, as a consequence, states with maximum 
hi- and tri-partite entanglement can be generated from propagation of unentangled photons. We 
generalize the results to n photons and show that entangled states with permutation symmetry 
can be generated from propagation of unentangled states. Generation of n-photon GHZ states is 
discussed as an example of a class of states with the desired symmetry. 

PACS numbers: 03.67.-a, 03.67.Bg, 42.50.-p, 42.50.Dv 



I. INTRODUCTION 

It is well-known that the classical coherence properties 
of an electromagnetic field vary due to propagation 
The most well-known example of this is the increase of 
the spatial coherence of the radiated field from an inco- 
herent source upon propagation 0, other examples 
are the Wolf effect, the variation of the spectrum of light 
under propagation 0, and the change of polariza- 
tion of light under propagation Q. It is therefore rea- 
sonable to pose the question: can propagation alter the 
quantum correlations of a field? In this paper, we study 
the change of entanglement on propagation and answer 
this question in the affirmative. A direct consequence 
of this result is post-selective generation of polarization- 
entangled photons through propagation of unentangled 
photons. For two photons this result has been pointed 
out by Lim and Beige Q and is similar in spirit to en- 
tanglement generation schemes in linear optics [l^l , 
where erasure of which-path information leads to cre- 
ation of entanglement. The scheme has also been im- 
plemented in reverse to generate entangled atoms by de- 
tection of photons [ll|j Il3| • Here, we extend these con- 
siderations to three photons and demonstrate that prop- 
agation and post-selective measurement can be used to 
create states with maximum genuine tri-partite entangle- 
ment [l^. Furthermore, we show that a generalization 
of this result leads to creation of n-photon Greenberger- 
Horne-Zeilinger (GHZ) states. 

Multi-photon entangled states have been generated 
for up to six photon by down-conversion and linear- 
optics 1141 and are of interest for optical quantum com- 
puting [1^, [l^. Moreover, many-particle entangled 
states are a resource in the one-way quantum computing 
paradigm [13, [3 which can be implemented advanta- 
geously in a linear optics setting [l^, [l^, HH . Interfer- 
ometric stability and beam-splitter alignment are major 
obstacles in creation of larger entangled states by lin- 
ear optics techniques. It would therefore be desirable 
to create multi-photon entanglement by simpler optical 



arrangements which may relax the requirement for beam- 
splitter alignment and stability. The schemes considered 
in this paper rely solely on free-space propagation and 
detection; our study offers insight into generation and 
manipulation of optical entanglement without the need 
for beam-splitters or non-linear optical elements. One 
major drawback of the proposed scheme is the exponen- 
tial scaling with the number of qubits due to the n-photon 
coincidence count. However, this deficiency can possibly 
be overcome via classical interference, and is currently 
under further investigation. 



II. TWO PHOTON CASE 

A. Generation of Entanglement 

As an example of a simple situation in which propaga- 
tion can change the quantum coherence of light, consider 
the situation of Fig[T] Suppose that one photon emerges 
from pinhole 1 in the polarization state a\H)-\-h\V), while 
a distance d from pinhole 1 a second photon is radiated 
from pinhole 2 in the state c\H) + d\V). Their combined 
state is therefore {a\H) + b\V)) ® {c\H) + d\V)). Two 
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FIG. 1: Two photons are emitted from two pinholes 1 and 2 
separated by a distance d. Two detectors 3 and 4 are posi- 
tioned in the far-field of the photons. There are two possible 
interfering processes by which the photons can reach the de- 
tectors. The interference between these paths leads to gener- 
ation of entanglement. 

point-like detectors 3 and 4 are positioned in the far-field 
of the pinholes. The photons reach the detectors after 
a time r/c, where their joint state (polarization) is mea- 
sured. We are only interested in the events in which both 
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detectors register a photon. For this geometry, the state 
of the photons at the detectors is pure and is given by 

1 pik{Ri3+R2i) 

HH)+b\V))^{c\H)+d\V)) 



-ic\H)+d\V))^ia\H)+b\V))} (1) 
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where k = lo/c is the wavenumber of the photons, Rij 
is the distance between source i and detector j and N 
is the normahzation constant. There are two interfering 
processes in which both detectors could register a single 
count: photon 1 landing on detector 3 and photon 2 on 
detector 4; photon 1 arriving at detector 4 and photon 2 
at detector 3. The two terms in the sum may be inter- 
preted as the two possible paths taken by the photons. 
In the far field of the sources where R13R24 — R23R14, 
the state vector simplifies to 



Ix) = ^{e-^''{a\H)+b\V))^{c\H)+d\V)) + 



e"^{c\H) +d\V))(^{a\H) +b\V))} 



(2) 



where (p = |(-Ri4-I--R23 — ^13 — -R24) and we have assumed 
that the radial terms in the denominator of Eq. ([1]) are 
all of the same order in the far field and can be absorbed 
in the normalization A''. The entanglement of this state 
can be quantified by concurrence [24I and is given by 
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2 cos2 (p(|ac|2 + + {\ad\^ + |c6|2 + 2\abcd\ cos ip) ' 

(3) 

As an example, if the photons begin in the state \HV), 
i.e. a = d = 1 and 6 = c = 0, one finds that the state at 
the detectors is -^{\HV) + \ VH)) and has entanglement 
of unity. 

This apparently counterintuitive result occurs because 
both photons are radiated into a large solid angle: 
whether a photon landing on a detector originated at 
1 or 2 is therefore unknown prior to measuring its polar- 
ization. We then post-select only those events in which 
both detectors register a count, projecting the detected 
state into a maximally entangled state. 



B. General Case 

The above analysis can be generalized for an initial 
state of the form a\HH) + P\HV) + -i\VH) + 6\VV), 
i.e. an arbitrary pure state of the two incident photons, 
straightforwardly. The concurrence is now found to be 
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1/?^ + 7^ + 2/37 cos 2lp - Aa5 cos^ ip\ 
1 + cos (2(^)(1 - |/3|2 - |7|2 + 2Re[(3r*]) 



(4) 



As a first example, the entanglement of the initial state 
-^{\HV) -\- \VH)) is invariant under propagation and 

the concurrence remains constant unity. The state 



^(|ffV) +i\VH)), however, generates a state with con- 
currence of |cos2(p|. Figure [2] shows a plot of concur- 
rence versus the phase Lp for the initial states \HV) and 
'^{\HV) +i\VH)). The entanglement of the initial state 

'^{\HV) + i\VH)) is destroyed and revived as the path 
difference between the two possible paths for the photons 
to reach the detectors, is varied. 




FIG. 2: Far field entanglement versus the phase ip, for 
the initial states -^{\HV) 4- i\VH)) (solid green line) and 
|//'F)(dashed blue line). Concurrence of the state generated 
from the initial state -^{{HV) +i\V H)) takes the simple form 
C{p) = \ cos{2tp)\ in the far field. (Color online) 



C. Beam-Splitter Analogy 

Our scheme is reminiscent of standard entanglement 
generation schemes in linear optics where two photons 
are incident on a beam-splitter and the events in which 
the photons are separated into two different ports are 
post-selected (Fig. We use the notation mj)AB = 




BS 



FIG. 3: The state emerging from the beam-splitter is entan- 
gled provided that the photons are separated into two ports 
A and B. 

\ni)A\n^j)B to represent a two-photon state in which 
n(m) j-polarized (j-polarized) photons are in the spa- 
tial mode A(B), and each of the subscripts i and j can 
take two possible values: H or V. For the incident state 
\^Ht^v)ab the state emerging from the beam-splitter 



3 



takes the form [22 



by (26 



N 2 



= -^{\lH,0)AB+i\0,lH)AB)'E) 

■^i-i\lv,0)AB + \0,lv)AB) (5) 



which is a product state and therefore not entangled. 
However, if we select the events in which two pho- 
tons separate into the two output ports A and B, 
the resultant state is the maximally entangled state 
-^{\^hAv)ab + \^vAh)ab)- Similarly, free space 
propagation acts like a beam-splitter, combining the 
hotons. The lack of the Hong-Ou- 
in free space propagation, however, 
means that there is a subtle difference between the two 
schemes. To see this, consider a general state of the form 

a\lH, '^h)ab+P\Ih, ly)As+7|lv, '^h)ab+S\1v, '^v)ab 
incident on a beam- splitter. Post-selecting the events 
at which each detectors registers a count, we find the 
detected state to be ^ , 1v)ab + 1//)ab), pro- 



state of the two 
Mandel effect d 



vided that at least one of the amplitudes /3 or 7 is non- 
zero. If both P and 7 are zero, photon bunching prevents 
the photons from arriving at different ports. This is in 
sharp contrast to the case considered in the previous sec- 
tion where the entanglement of the detected state has a 
strong dependency on the amplitudes of the initial state. 

Similar post-selective schemes have been utilized in 
linear optics quantum computing [l^ . [23 | to gen- 
erate entangled photons and perform logical operations 
or to induce effective non-linearities [2^. Since all such 
schemes rely on beam-splitters to erase the which-path 
information, their scalability is severely limited by in- 
terferometric stability; it would, therefore, be desirable 
to create multi-photon entangled states by simpler ar- 
rangements that do not rely on beam-splitters and/or 
Hong-Ou-Mandel effect. 



III. TWO PHOTON CASE: EMISSION BY TWO 
ATOMS 



Having established the role of interference between 
paths in entanglement creation/destruction, let us now 
consider a concrete example where the photons originate 
from spontaneous emission of two atoms. Apart from a 
means of implementing the scheme experimentally, this 
calculation places the heuristic argument presented above 
on firmer ground, and allows us to compute the entan- 
glement at the intermediate points from the near-field to 
the far-field of the source. For a system consisting of two 
three-level atoms, interacting with a quantized field, the 
Hamiltonian in the rotating wave approximation is given 



i X 



k.AOkAS^j + C.C 



(6) 



where J2i indicates summation over the atoms, J^k is the 
vector sum over the spatial field modes, J2\ i^ the sum- 
mation over the two orthogonal polarizations, s\i {s\^) 
is the atomic lowering (raising) operator acting on the 
i*^ atom and corresponding to a transition with a A po- 
larized photon, OkA (ajkA) field annihilation (cre- 
ation) operator and 5i,k,A is the couphng constant. The 
two excited states are assumed to be degenerate with luq 
the atomic transition frequency. We choose to work in a 
basis where the atomic operators correspond to linearly 
polarized photons. Solving the Heisenberg equation of 
motion, ihsxiit) — [sxi{t), H], we arrive at the follow- 
ing differential equation for the slowly varying amplitude 



r,iUl0t 



dt 2 ^ y "J w 

i 



(7) 



where the matrix My describes the interaction between 



two atoms separated by d, 



M, 







(8) 



The solution to the differential equation is given by 
s„,(t) = e-^'^«*^O^^^Cn^^^e'^"''*/'e-''(i+'^-)*/2s„,(0)(9) 

where Cn*'^'' is the n'^ eigenvector of the interaction ma- 
trix A/, with eigenvalue a-p ~\- iPp. In physical terms ap 
and (3p are the frequency shift and the decay shift due 
to super-radiant effects which are small if the atoms are 
more than one wavelength of radiation apart (27j and 
may also be neglected. As a consequence of the orthog- 
onality of the eigenvectors [2^, the Heisenberg equa- 
tion for the atomic operators is reduced to Soii{t) — 
■Sai(0)e~*"°*e~'>'*/^ for a = 1 or 2, which is equivalent 
to the solution of a classical oscillating dipole. The elec- 
tric field at detector j in the semi-classical approximation 
is given by the formula |29| 



+ 3(nj[n.s,(t; .) - s;(t; .)]) 



'nj X [nj X s,(i, .)] 



ik 



(10) 



where j = {3,4}, Rij is the distance between atom I 
and detector j, is the retarded time, = t — Rij/c, 
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§1 = {su, §21,0), p is the amplitude of the dipole moment amphtudes of the detected state are therefore given by 



of the atom and rij is the unit vector in the direction of 
the observer. We have chosen the coordinates such that 
Sii{s2i) is the component of the dipole moment operator 
along the x-axis (y-axis) (Fig. [4]). We use \x), \y) and 
\g) to denote the eigenstates of the atomic operators, we 
can therefore write Siisi2|a;a;) = \gg), S2iSi2|ya;) = \gg) 
and so on. The radiated field may be simplified for an 
observer in the far-zone; however, since we are interested 
in computing the propagation of entanglement, all the 
terms in Eq. pil)l will be kept. We compute the state vec- 
tor of the post-selected two-photon state reaching the de- 
tectors for symmetric detection angles. The components 
of the electric field at detector j decomposed along the az- 
imuthal unit vectors eg = (cos cos 6*, sin cos 6*, — sin 6*) 
and = (— sin (j>, cos 4>, 0) are given by 



Sxi{t) sin (/)) 



1=1 

1 ik 



Rij Ru 



(11) 



(+), 



(13) 



where rj and k can he 9 or (j) and \gg) is the ground state 
of the atoms. 

As an example for the initial atomic states {ipi) = l^u), 
\ip2) = -^{\xy) + \yx)) and IV-s) = -j^{\xy)+i\yx)), where 
\x) and \y) represent the state of an atom in terms of the 
direction of its dipole moment, the detected two-photon 
states are respectively given by 



lxi) = ^(|i?y) + e- 



-2ikd sin 6 



\VH)) 



1X3; = 



1X2) = ^^\HV) - \VH)) 



~2ikd £ 



\VH). 



Sxl{t) COS 0008 6*) 



' 1=1 
P 1 



{syi{t) COS 9sm(j) + 
ik 



R13 Rij 



Ru 



g^fcK,, (12) 



where s^i = Si/ and Syi = §21 and the retarded 
dipole moment operator has been written as 
s{t — Rij/c) — s{t)e'^^^^' . We can simphfy the 





FIG. 4: Two single photon sources 1 and 2 and two detectors 3 
and 4 are positioned in the xz-plane as shown. The detectors 
are positioned symmetrically about the z-axis. The photon 
polarization is measured along the unit vectors ee and 60. 

calculations by remembering that all detected two- 
photon states are pure. A general (unnormalized) 
state reaching the detector is therefore of the form 
Ix) = Ce^\HV) + C^g\VH) + Cee\HH) + C^^\VV) 
where \H) and \V) represent the two linear po- 
larizations in the detector basis, along the unit 
vectors eg and respectively. From the the- 
ory of photo-detection [33| we know that jCe^p — 

{^\Ei;^ (rns)^;) irn^)Eit^ {rn^)E%^ (rng) |^) where 
\'ip) is the initial state of the atoms and, since t = r/c, 
we have omitted the time dependency of the fields. The 
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FIG. 5; Concurrence versus the detector position r, in units of 
\/2n for 9 = 45°. The initial atomic states are |a:i/) (dashed 
blue line) and -^{Ixy) + i\yx)){so\\d green line), a) kd — 
3.45 (far-field minimum) and h) kd = 7 (far-field maxi- 
mum). (Color online) 

Ixi) and 1x2) have a concurrence of unity and jxa) has 
concurrence of | cos {2lp) \ , in agreement with the heuris- 
tic treatment of the previous section. In fact, if the ini- 
tial state of the atoms is an arbitrary state of the form 
a\xx) + (3\xy) + j\yx) + 5\yy), one arrives at Eq. ([4]) for 
concurrence in the far-field. 

We now compute the entanglement at all intermediate 
points as the detectors are moved from the near-field to 
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the far-field. For a detection angle of 45°, the results 
are presented in Fig. ([5]). The behavior of the state \xy) 
is particularly intuitive to understand; at r = the de- 
tectors are at the origin and concurrence is unity due to 
complete mixing of photons. A "near-zone" minimum oc- 
curs at r ~ I where a detector is positioned close to each 
atom and the probability of a photon captured by the 
farther detector is negligible. Concurrence then recov- 
ers its far-field value as the detectors are moved further 
apart. The state \xy) + i\yx) also has unit concurrence 
at r = 0, but no subsequent "near-zone" minimum. This 
is expected since the initial state was maximally entan- 
gled and the entanglement is directly transferred from the 
atoms to the photons if a detector is placed next to each 
atom. For r < d/2 a series of interference fringes occur 
as the detectors are moved apart. Concurrence recovers 
the far-zone value for r >> d. 
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FIG. 6: Concurrence for symmetric detector positions for the 
initial atomic states \xy) {dashed blue line) and -^(|a:j/) -I- 
i\yx)){so\id green line). Both detectors make an angle d with 
the z-axis, where — 7r/2 < d < n/2. a) Concurrence in the 
near field for kd = 4.45 and r — 2. b) Concurrence in the far 
field for kd = 4.45 and r = 1000. (Color online) 

Finally to see the spacial variation of entanglement, 
we move the detectors in the xz-plane, keeping them 
symmetric at all times, and plot the variation of con- 
currence as a function of the detection angle 6 (Fig. [4]). 
Figure [6^. shows the entanglement in the near-field. Fig- 
ure [Hb- is the corresponding plot in the far-field where 
we have chosen the atomic separation such that the state 
-^{\xy) + i\yx)) would have maximum concurrence for a 

6 — 45° detection angle. The state \xy) generates max- 
imally entangled photons for all symmetric detector ori- 
entations. Similar fringe patterns are observed at both 



extremes for the state -^{\xy} + i\xy)), confirming the 
previous observations that altering the path difference 
between the photons can create maxima and minima of 
entanglement. Our results are consistent with the find- 
ings of Lim and Beige [7|] who study the spacial variation 
of entanglement of formation for two dipole sources, but 
do not discuss the variation of entanglement at the inter- 
mediate points. 



IV. N-PHOTON GHZ-STATES 

A. General Result 

The above observations give rise to the question: can 
propagation and post-selection be used to create multi- 
photon entangled states? In this section, we first state 
a general symmetry property of any n-photon state that 
can be generated via propagation and post-selection in 
our chosen geometry. We then show how the initial con- 
ditions can be tailored to create n-photon GHZ states up 
to local unitaries. Finally we consider a three-photon 
GHZ state as a specific example. We assume an ar- 
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FIG. 7: We consider an array of single photon emitters pre- 
pared in the state I"!/)) = X^i^ i i'ij^...^^ |ii...i„). Polarization- 
resolving detectors are positioned in the far-field of the atoms 
with the same spacing as the atomic lattice. 

rangement consisting of a one-dimensional array of sin- 
gle photon emitters, separated by a distance d, with a 
single photon detector directly above each emitter a dis- 
tance z away (Fig. [7]). We consider an initial state of 
the form lip) = Z]ii...i„ Ai...z„ |«i---«n)- After propaga- 
tion and post-selection, this produces a final state of the 
form Ix) = i C'ii...i„ Ki---*n) where, each of the sub- 
scripts can take two possible values, H = or V = 1. 
We demonstrate that states with permutation symmetry 
can be generated from propagation of a suitable initial 
state, provided that the detectors are far enough from 
the source. Local unitaries may need to be applied on 
the photons before detection to create the desired entan- 
gled state. A state |x) = Ci^,,,i„\ii---in) possess 
permutation symmetry if all amplitudes Ci-^,,,i^ with dif- 
ferent permutations of the subscripts are equal. The far- 
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field condition is the key behind this result and assumes 
a first order approximation for the distance between each 
emitter and each detector; i.e. the phase difference be- 

2 .2 

tween different paths is neglected. This demands ^ 
1 and becomes more difficult to satisfy for large n. The 

n-photon GHZ state \GHZ),, = {\H)®" - 1^)®"), 
is an example of an state with the desired permutation 
symmetry; it can thus be created by simple spatial prop- 
agation, if the amplitudes of the initial state are chosen 
carefully. We demonstrate in Appendix A that propaga- 
tion of the separable state 

IV^) = (g) (^sin (^^) \H) + cos (^^) IF)) (14) 

generates the post-selected entangled state |x) up to 
global phases in the far-field such that 

(H5t)^"|x) = -i=(|i?)®"-|y)«") (15) 

where H is the Hadamard gate and S is the phase gate 
such that S\H) \H), S\V) = i\V), n\H) = -^{\H) + 

\V)) and n\V) = ■^{\H) - \V)). Eq. ^ corresponds 
to an initial state of the atoms such that the polarization 
for the n-th emitted photon simply forms an angle of 
with the horizontal. 

Assuming emission into a 47r solid angle and a total 
detection solid angle of fi, the probability of an n photon 
coincidence is Pn = nlrj'^ for rj — ijqTia, where ijq is the 
quantum efficiency of each detector and rj^ = Q/ (Ann) 
(assuming that the area of each detector scales as i). 
For large n, using Sterling's formula for large factori- 
als, one can show that the n-photon coincidence scales 
as Pn ~ -\/27m (rygJ7/(47re))"; in other words, it scales 
exponentially with the size of the array. The main draw- 
back of the scheme, is therefore, the low probability of 
registering an n-photon coincidence. One possibility of 
overcoming this obstacle is to use classical interference 
to maximize the probability of detection at the desired 
detector locations and will be investigated in a subse- 
quent publication. 

B. Three Photon GHZ state 

As a concrete example and to elucidate Eq. and 
Eq. p^ . consider a general three-photon state of the 
form |-0) = X^rj/te -^'/keI^kc) which generates the state 
Ix) = SrjKc C''7K«l'7«e) in the far-field. By generalizing 
Eq. (|13p the amplitudes of the detected state can be ex- 
pressed as 

C,,,, = {ggg\E^+\rn^)Eif{rns)El^t\rnem. (16) 

This equation can be simplified for an observer in the 
far-field. Inserting Eq. ^ and Eq. ^ into Eq. HH) 



for = 0, we arrive at the far-limit of Eq. (fT6|) 

Ctik^ = {ggg\ ^ ^ Dap^ |a/37). (17) 

Imn a/37 

For the arrangement considered Rji = z for all i and j in 
the far-field and therefore all radial terms have dropped 
out. The amplitude C^^^ is therefore the sum of six 
terms: all cyclic and anti-cylic permutations of 77, k and 

Cjy/tf = -D^K? + -Dfc^j) + D^r]K. + Djj^K. + Dnn^ + -Dffcj) (18) 

This immediately proves, for example, that the state 
\HHV) generates a W state in the far- field. To see 
how GHZ states are created, one must demonstrate that 
propagation of the initial state \tp) = {^f=i(sin(^)|i/) + 
cos(^)|T^)) generates the state |x) = Er^Ke C',,„e|f?Ke) = 
\{\HHV) + \HVH) + \VHH) - \VVV)) in the far-field. 
To prove this it is sufficient to show that a) Cqoi = 
Coio — C'loo = —Cm, b) Cuo = Cioi = Con = . 
Both these criteria can readily be verified from Eq. (|18p . 

Three qubit GHZ states have genuine tri-partite entan- 
glement [13| and show maximum violation of three-qubit 
Bell inequality [3l| . Three-tangle is a measure of genuine 
tri-partite entanglement and is defined to be (l3| 

— ^1(23) ~ ^12 ~ ^13 (19) 

where Ctj is the concurrence between qubits i and j and 
Ci(23) measures the entanglement between qubit 1 and 
the joint state of qubits 2 and 3. Three-tangle is bounded 
between and 1 and states equivalent to GHZ states up 
to local unitaries are characterized by a three-tangle of 
unity. Figure [5] is a plot of three-tangle versus distance 
for the initial state |^) = {8)r=i(sin(^)|i?) -l-cos(^)|l/)) 
for two different initial atomic spacings. In the far-field 
all phase information is washed out and three-tangle ap- 
proaches unity independent of the initial separation as 
expected. 




z 



FIG. 8: Three-tangle versus the perpendicular detector po- 
sition z in units of \/2n for the initial state described by 
Eq. JTHl (n=3) for d = A (solid red line) and d = 2A (dashed 
blue line). 
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V. CONCLUSION 

In conclusion, we have demonstrated that entangle- 
ment of photons emitted in a large solid angle can change 
on propagation. This arises because in the far-zone 
all information about the origin of the photons is lost 
and this leads to quantum mechanical interference be- 
tween all possible paths to the detectors. We use con- 
currence and three-tangle as measures of two and three 
qubit entanglement and verify that both these quan- 
tities vary smoothly from near-field to far-field. We 
demonstrate that the propagation of the state = 
<8>r=i(sin(^)l-ff) +cos(^)|T/)) generates n-photon GHZ 
states post-selectively. Our results appear to suggest that 
the chosen geometry is suitable for generation of states 
with permutation symmetry; both n-photon W and GHZ 
states fall into this category. 

The scheme may be realized experimentally using 
quantum dots or ions prepared in arbitrary initial states. 
The main drawback of the scheme is that the entan- 
gled photons are only accessible via post-selective or non- 
demolition measurements, which makes the scheme, in its 
current form, of limited use in practical generation of en- 
tanglement. However, generation of effective interactions 
between photons without the need for beam-splitters has 
obvious attractions in the design of linear optics quan- 
tum processors. One possibility of overcoming the post- 
selectivity criterion is using classical interference to max- 
imize the probability of photon detection at the desired 
detector locations. This possibility is currently under in- 
vestigation. 
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APPENDIX A: n-photon GHZ-states 

Here we present the proof that for a particular choice 
of the initial state, the state reaching the detector is 
an n-photon GHZ state up to local unitary rotations. 
We adopt the binary notation to denote the polarization 

states, 

|0) ^ \H) (20) 

|1) ^ \V). (21) 

Theorem: For the specific geometry shown in Fig[7| 
propagation of the state 




post-selectively generates the state |x) up to global 
phases in the far-field such that 

(H5t)®"|x>^-^(|0)«"-|l)«") (23) 

2 ,2 

where the far- field is defined by the condition ^^p- ^ 
TC is the Hadamard gate and S is the phase gate such 
that S\0) = |0) and S\l) = n\0) = ^(|0) + |1)), 

Till) = ^(|0) - 

Proof: Let us assume initially that n is odd. The initial 
state and the final state |x) can be rewritten as\tp) = 
Eii...»„ Ai...»„Ni-.in) and \x) = Eii...»„ C'ii...,„|ii...i„). 
By generalizing Eq. we have 

ni perm. 

where the summation is carried out over all n\ possible 
permutations of the subscripts. The detected state pos- 
sesses permutation symmetry; meaning that all ampli- 
tudes with different permutations of the same subscripts 
are equal. We use the notation C„^m to denote the ampli- 
tude of an n-photon state in which m are in the state 1. 
We must show that the detected state with no rotation 
is of the form 

{snr" ^ (|oo...o)-|ii...i)) = 

E (25) 

il...in 

where the amplitudes are symmetric with respect to per- 
mutations and can be expressed as 

^ f m even , , 

(— 1) 2 ra odd. 

The strategy is to prove that Cn,m = Bn,m- In order to 
do so we must show that: a)Cn,m — for all even m, b) 

C m — l 

g"'™ = (— 1)^~ for all odd m with n > m. 

Condition a. The amplitudes Di-^ ^ i^ can be written 
as Di^ i^^ — H^l" /('^•ij) where the function f{x.ij) is 
defined to be 

. . / sin(a::) if ij = 

= I cos(4 if *J = 1 (27) 

We therefore have = E„! pe™. 0^=1 / • 

Since i„ = 1 we can exclude the last photon and write 

c^^-^- E n/(^-^.) (28) 

(n — 1)! perm, j — 1 

which vanishes if the product contains an odd number of 
cosines, or an even number of Is (including the in pho- 
ton). To see this consider the amplitude C„^2- Expanding 
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Eq. p8|) we arrive at 

Cn.2=l Vsin(^)cos(-) n 



. n n 

n—1 . . n—1 . \ 

5: sin (^) cos (^) n Bin^(^))(n-2)!2! 



Condition b therefore demands: 



E 

il>i2>---i„ 



cot^ ( — )... cot^ (:^) ^ ^^\, (32) 
rt 71 m\(n — my. 



which vanishes since sin(— ) = sin(22_^7r) and cos(— ) = 
— cos(^2^7r). This proves condition a. 

Condition b. We now write expressions for C„ i and 
Cn^rm excluding all vanishing terms in Eq. (j28p we obtain 



C„4 = (n-1)! ji^sh 



(29) 



This equation can be proven by recalling that tr = 
cot^ (2pTt) ^ ~ li2, ...p are the distinct roots of the 
pth degree polynomial [3^ 



and 



n — 1 

C„,„ = (-l)"'(n - m)!m! ^ cos^ 



COS 



(^) n 



2l>«2>---« 

kn 



sm 



k^i^ .a—l...m' 



n 



(30) 



where m' The ratio of Cn,i and C„_m is therefore 

given by 

n-l 

(1 2 / ■ 



cot^ 



n 



ii >i2--->i„i/ 

2 /'*m'7r\ — m)!m! 



cot" 



n y (n — 1) 



(31) 



Q(t) 



fe=0 ^ ^ 



(33) 



where 



fe!(n-fc)! 



By applying Viete's formula (3E 



to Eq. p3|) and making the substitution p = 22__L and 



m 



m — 1 
2 



we arrive at Eq. 



This proves condition b. 



If n were assumed to be even from the outset, the up- 
per limit of the sum in Eq. ((32)l would have been § — 1- 
The proof for even n would have then demanded the sub- 
stitution p = f — 1 in Eq. [551 
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